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Within an effective Lagrangian model, we present calculations for cross 
sections of the dilepton production in proton-proton and proton-neutron colli- 
sions at laboratory kinetic energies in 1-5 GeV range. Production amplitudes 
include contributions from the nucleon-nucleon bremsstrahlung as well as from 
the mechanism of excitation, propagation, and radiative decay of A(1232) and 
N*(1520) intermediate baryonic resonances. It is found that the delta isobar 
terms dominate the cross sections in the entire considered beam energy range. 
Our calculations are able to explain the data of the DLS collaboration on the 
dilepton production in proton-proton collisions for beam energies below 1.3 
GeV. However, for incident energies higher than this the inclusion of contri- 
butions from other dilepton sources like Dalitz decay of 7r° and r\ mesons, and 
direct decay of p° and to mesons is necessary to describe the data. 
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I. INTRODUCTION 



Dileptons observed in the nucleus-nucleus collisions travel relatively unscathed from the 
production point to the detector. Therefore, they are expected to provide clear information 
about the early dense and hot stage of heavy ion collisions [1-3], which is in contrast to 
the hadronic probes which often suffer from strong final state interactions and the infor- 
mation about the collision history carried by them may be lost due to the rescattering in 
the expansion phase [4]. One of the phenomena predicted [5-13] at higher nuclear matter 
densities is the restoration of chiral symmetry which is manifested in the modification of 
masses of the vector mesons as a function of the nuclear matter density. Consequences 
of this effect can be observed in the dilepton (e + e~) spectra measured in nucleus-nucleus 
collisions [14]. Enhancements (above known sources) observed in the measured [15,16] soft 
lepton pair production cross sections in heavy ion collisions at the SPS energies, have been 
attributed [17-21] to in-medium modifications of the vector meson properties rather than to 
reflection of a new state of hadronic matter. Although, such a scenario has not been fully 
successful [22-24] in explaining the dilepton spectra measured in these reactions at much 
lower beam energies (a few GeV/nucleon) [25] where the temperature and density regime is 
different and one does not expect a phase transition. 

In this context, the investigation of the dilepton production in elementary nucleon- 
nucleon (NN) collisions is of interest because the corresponding cross sections enter in 
the transport model calculations of the e + e~ spectra in the heavy ion collisions. Therefore, 
a quantitative understanding of this process is a natural prerequisite to an unequivocal de- 
termination of in-medium effects mentioned above. The study of this process is of interest in 
its own right as it is expected to provide deeper insight into the hadronic structure and the 
photon-hadron interactions. This investigation may give fundamental information on the 
electromagnetic form factor of the nucleon in the time-like region around the vector-meson 
masses which is otherwise hard to access(see, e.g., [4]). 

It has been known for some time that intermediate baryonic resonances play an important 
role in the dilepton production in the NN collisions. Dalitz decay of the A isobar has 
been shown to be a strong dilepton production channel [26-28]. The importance of the 
baryonic resonance N*(1520) has been emphasized for this process in Refs. [29,30] where 
it is pointed out that the subthreshold p production (and its subsequent decay) via this 
resonance makes important contributions to the dilepton spectrum observed in the proton- 
proton (pp) collisions. Therefore, the investigation of the dilepton production is expected to 
provide a useful tool to probe the parameters (e.g., coupling constants, form factors etc.) 
of the nucleon-resonance-photon vertices and masses and widths of the relevant nucleon 
resonances. 

The aim of this paper is to investigate the dilepton production in nucleon-nucleon col- 
lisions in the beam energy range of 1-5 GeV within an effective Lagrangian model (ELM) 
which has been successfully applied to the description of pion and associated kaon produc- 
tion in pp collisions [31-33]. Initial interaction between two incoming nucleons is modeled 
by an effective Lagrangian which is based on the exchange of the 7r, p, uo and o mesons. 
The coupling constants at the nucleon-nucleon-meson vertices are determined by directly 
fitting the T matrices of the nucleon-nucleon (NN) scattering in the relevant energy region. 
The effective Lagrangian uses the pseudovector (PV) coupling for the nucleon-nucleon-pion 
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vertex and thus incorporates the low energy theorems of current algebra and the hypothesis 
of PCAC. The e + e~ production proceeds via excitation, propagation and radiative decay of 
A(1232) and A* (1520) baryonic resonance states. Also included are the nucleon intermedi- 
ate states (which gives rise to the AW bremsstrahlung contribution). The interference terms 
between various amplitudes are taken into account. The gauge invariance at the electro- 
magnetic vertices is preserved in our calculations. Our model is similar in spirit to those of 
Refs. [28,34,35]. However, in Ref. [35] no resonance contribution was considered whereas in 
Ref. [28,34] resonance contributions were limited to the A isobar only. The latter, though, 
did include the interference between the nucleon and A terms. In calculations presented in 
Refs. [29,30,36], AW bremsstrahlung contributions were ignored. We would like to stress 
that in Refs. [29,30] constant matrix elements have been used for various processes and the 
total dilepton production cross sections have been calculated by adding the corresponding 
cross sections and not the amplitudes, so part of the motivation of the present study is the 
investigation of the so far neglected quantum mechanical effects. 

We investigate the role of baryonic resonances in the invariant mass spectrum of the 
dilepton produced in proton-proton and proton-neutron collisions at various beam energies 
in the 1-5 GeV range. To this end, we present the first field theoretic calculation of the 
dilepton production in AW collisions where excitation, propagation and radiative decay of 
the N*(1520) baryonic resonance are fully accounted for. We also compare our calculations 
with the published data [37] on the dilepton production in elementary proton-proton colli- 
sions, by the Dilepton Spectrometer (DLS) collaboration. In order to describe these data 
the contributions from other dilepton sources (n° and rj Dalitz decay and direct decay of p° 
and uj mesons) have also been considered. 

The remainder of this paper is organized in the following way. Section II contains details 
of our theoretical approach. Section III comprises the results of our analysis and their 
discussions. The summary and conclusions of our work is presented in Sec. IV. 

II. FORMALISM 

A representative of the lowest order Feynman diagrams contributing to the dilepton 
production as considered by us, is shown in Fig. 1. The intermediate nucleon or resonances 
can radiate a virtual photon which decays into a dilepton (Figs, la and lb). There are 
also their exchange counterparts. In addition, there are diagrams of these types where the 
virtual photon is emitted from the nucleon line on the right side. The internal meson lines 
can also lead to dilepton emission (see Fig lc). Momenta of various particles are indicated in 
Fig. la. q, pi, and k are four momenta of the exchanged meson, the intermediate resonance 
(or nucleon) and the photon, respectively. To evaluate various amplitudes, we have used 
the effective Lagrangians for the nucleon-nucleon-meson, resonance-nucleon-meson, nucleon- 
nucleon-photon and resonance-nucleon-photon vertices. These are discussed in the following 
subsections. 



3 




<-- 

q 



p >/\ p. a 



(a) 





(b) 



(c) 



FIG. 1. A representative of Feynman diagrams for emission of dilepton in nucleon-nucleon colli- 
sion as considered in this work, (a) denotes emission after NN collisions, (b) before NN collision 
and (c) during NN collision. The box represents any of an off-shell nucleon, a A isobar or a N* 
resonance. 



A. Nucleon-nucleon-meson vertex 

As done before in the investigation of pp — > ppir , pp — > pmr + [31] and pp — > pAK^ 
reactions, the parameters for these vertices are determined by fitting the NN elastic scatter- 
ing T matrix with an effective NN interaction based on the n, p, u, and a meson exchanges. 
The effective meson- AW Lagrangians are (see, e.g., [38,39]) 

CNN* = -f^7V757^T • (d^)* N . (1) 

IrriN 

Cnnp = -9nn p ^n (lp + 2~~ < V 9 '') r ' P ^ N - W 

£nNw = -gNNco^N (ip + ^~ (T l»' a ^j U ^N- (3) 
£>NNa = gNNa^NO-^N- (4) 

In Eqs. (l)-(4), we have used the notations and conventions of Bjorken and Drell [40] and 
definitions of various terms are the same as those given there. In Eq. (1) tjin denotes the 
nucleon mass. It should be noted that we use a PV coupling for the NNtt vertex. Since 
these Lagrangians are used to directly model the AW T matrix, we have also included a 
nucleon-nucleon-axial-vector-isovector vertex, with the effective Lagrangian given by 
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Cnna = gNNA^lnl^ ■ A>", (5) 

where A represents the axial- vector meson field. This term is introduced because if the mass 
of the axial meson A is chosen to be very large (^> tun) [28] and Qnna is defined as 

g NNA = ^—j , (6) 

it cures the unphysical behavior in the angular distribution of NN scattering caused by the 
contact term in the one-pion exchange amplitude. It should be mentioned here that A is 
different from the ai(1260) meson resonance to be discussed later on. The role of the A 
vertex is to explicitly subtract out the contact term of the one-pion exchange part of the 
NN interaction. Similar term in the coordinate space potential is effectively switched off by 
the repulsive hard core. 

At each interaction vertex, the following form factor is introduced 

Fr=(%^),i=ww (7) 

where qi and m 8 are the four momentum and mass of the ith exchanged meson and \ is 
the corresponding cut-off parameter. The latter governs the range of suppression of the 
contributions of high momenta which is done via the form factor. Since NN scattering cross 
sections decrease gradually with the beam energy (beyond certain value), and since we fit 
the elastic T matrix directly, the coupling constants are expected to be energy dependent. 
Therefore, while keeping the cut-offs \ [in Eq. (7)] energy independent, we take energy 
dependent meson-nucleon coupling constants of the following form 

g(V$) = goexp(-£y/s), (8) 

here s is the square of the total CM energy. The parameters go, £, and A were determined 
by fitting to the T matrix to the relevant proton-proton and proton-neutron scattering data 
at the beam energies in the range of 800 MeV to 4.00 GeV [28]. This procedure also fixes 
the sign of the effective Lagrangians [Eqs. (l)-(5)]. The values of various parameters are 
shown in Table I [the signs of all the coupling constants (g) are positive]. In this table 
the parameters of the A exchange vertex are not explicitly shown as they are related to 
those of the pion via Eq. (6). The mass of the A meson is taken to be 18 GeV. The same 
parameters were also used to describe the initial NN interaction in the calculations reported 
in Refs. [28,31,33]. This ensures that the elastic NN elastic scattering channel remains the 

TABLE I. Coupling constants for the iViV-meson vertices used in the calculations 



Meson 


g 2 /4n 


I 


A 


mass 








(GeV) 


(GeV) 


7T 


12.562 


0.1133 


1.005 


0.138 


a 


2.340 


0.1070 


1.952 


0.550 


to 


46.035 


0.0985 


0.984 


0.783 


P 


0.317 


0.1800 


1.607 


0.770 


k p = 6.033, = 0.0 
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same in the description of various inelastic processes within this model, as it should be. 

The exchange of a\ meson (mass = 1260 MeV), which is the chiral partner of the p meson 
has also been considered by some authors [41] in the description of the NN interaction. 
However, its contribution there is masked to a considerable extent by the strong short- 
ranged repulsion originating from u exchange. Moreover, as can be seen from Table I, the 
values of the cut-off parameters appearing in the form factors [Eq. (7)] are in the region of 
1.0 - 2.0 GeV. Therefore, it does not seem meaningful to include exchange of mesons with 
masses lying in the region where modifications due to the extended structure of hadrons 
are applied. Nevertheless, a\ meson (if included) can contribute to the dilepton spectrum 
through diagrams like the one shown in Fig. 1(c). However, the couplings of this meson to 
nucleons and resonances are largely unknown. In any case, contributions of such processes 
are expected to be of relevance only for the dilepton invariant masses > 1.5 GeV [42]. This 
region is excluded by the DLS data. Due to these reasons we have not included the exchange 
of ai meson in our calculations. 



B. Resonance- nucleon- meson vertex 



In addition to nucleonic intermediate states (NN bremsstrahlung) , we have considered 
in this work also the contributions from the A(1232) isobar and iV*(1520) nucleon resonance 
intermediate states. The latter is a spin-| negative parity resonance. While only ir and p 
mesons couple to A isobar, all of the four exchanged mesons, namely, n, p, u>, and a can 
couple to iV*(1520) resonance. 

For spin-| resonances, we use the following effective Lagrangians [43] 



Urn, = 9 -^^ R V^^ N + H.c. (9) 
Crn p = i^*J {d v ff - d»p v ) lv Y^ N + H.c. (10) 



m p 

c RNu , = i^^J (d^ - 7 ,r^ + H.c. (11) 

C RNa = 9 -^^ R T„{d»<j)* N + H.c. (12) 
m a 

Here, is the vector spinor for the spin-| particle. In Eqs. [(9) - (12)], the operator T n 

TABLE II. Coupling constants and cut-off parameters for the resonance-nucleon-meson vertices 
used in the calculations 



Resonance 


decay channel 


9 


cut-off parameter 
(GeV) 


A(1232) 


Nit 


2.13 


1.421 




Np 


7.14 


2.273 


iV*(1520) 


Ntt 


1.55 


0.800 




Np 


6.44 


0.800 




Nlu 


3.42 


0.800 




Na 


1.24 


0.800 
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is unity for even parity resonance and 75 for the odd parity one, whereas T is 75 for the even 
parity resonance and unity for the odd parity one. The meson fields in above equations need 
to be replaced by r • (where corresponds to n, p, u, or a meson fields) and T • (where 
T and r are the isospin transition operator and Pauli isospin matrices, respectively) for 
isospin-^ and isospin- 1 resonances, respectively. The values of the coupling constants g&N-w 
and QANp are shown in Table II; these are the same as those determined in Ref. [28] by fitting 
to experimental data on pp — > nA ++ reaction at kinetic energies in the range of 1-2 GeV. We 
have assumed that the off-shell dependence of the resonance-nucleon vertex is determined 
solely by multiplying the vertex constant by a form factor. Similar to Ref. [28,31], we have 
used the following form factor for the A vertices 



F NA = 



(Af f - mf 
(Af ) 2 - qf 



i = 7i,p, (13) 



with the values of the cut-off parameters as Af and Af as given in Table II which are exactly 
the same as those used earlier [28,31,33]. 

For the N*Nn and N*Np vertices the coupling constants have been determined from 
the observed branching ratios for the decay of the resonance to Nit and Np channels, re- 
spectively. In the later case the finite lifetime for the decay p — > nn has been taken into 
account by introducing an integration over the corresponding phase space. The details of 
this method are provided in Ref [33]. The coupling constants gN*Nw and Qn*n<j are deter- 
mined by vector meson dominance (VMD) hypothesis [44] and from the branching ratio of 
the decay of this resonance into the two pion channel [33], respectively. It should be noted 
that there is considerable uncertainty in the latter two coupling constants. However, the 
contributions of these terms to the total dilepton production amplitude are almost negligi- 
ble. As branching ratios determine only the square of the corresponding coupling constants, 
their signs remain uncertain in this method. Predictions from the independent calculations 
can, however, be used to constrain these signs [43,45,46]. Guided by the results of these 
studies, we have chosen the positive sign for the coupling constants for these vertices. The 
values of various coupling constants are given in Table II. 

As in Refs. [43,44], we have used the following form factor for the N*N vertices 



NN* 



(Af* 



(Af )4 + ( ? |_ m 2 )2 



,j = ir,p,u,a, (14) 



with the value of the cut-off parameter being 0.8 GeV in all the cases. It may however be 
mentioned that identical results will be obtained if one uses the form factor as given by 
Eq. (13) with different cut-off parameters. 

C. Nucleon-nucleon-photon vertex 

In the nucleonic bremsstrahlung process of dilepton production the intermediate nucleon 
is necessarily off-shell. The general form of the effective Lagrangian is given by 

C NNy = -e^rJ^ttA", (15) 
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where the half-off-shell nucleon-photon vertex function T u is [47-49] 

T NN, = _ (F s^ + ps^ + F ^ m)As . (16) 

s=± 

In Eq. (16), k 2 = (pf — pi) 2 with pf and being the initial and final nucleon four momenta. 
A± = (±^ + W)/W are the projection operators where W = (p 2 ) 1 / 2 and S M = icr uu k u /2m^. 
The form factors -^1,2,3 are functions of A; 2 , VT, and m^. The factor F 3 is not independent 
but is constrained by the Ward-Takahashi identity [50] (i.e. by the requirement of gauge 
invariance) 

F ' = i " + ^r/"' (17) 

where is the nucleon charge in units of | e |. 

In this work, we have adopted the procedure followed in Refs. [28,51] where one uses 
the on-shell form of the vertex function also for the off-shell momenta. This means, one 
assumes = F{ = F 1 and F 2 = F 2 = F 2 and F 3 + = F 3 ~ = 0. This vertex does 
not in general satisfy gauge invariance. However, this can be achieved (see also [52]) by 
multiplying the external photon emission vertices by the same form factors that multiply 
the hadronic vertices [Eq. (7)] and by multiplying the internal photon production diagrams 
by the following additional factor: 

2 22 2 

zp 1 1 meson HI , ''''meson V2 /i q\ 

F ^~ 1+ X 2 -ql + \ 2 -ql ' (18) 

where qi and q 2 are the four momentum transfers at the left and right vertices, respectively. 
In this way various electromagnetic form factors can be implemented for the hadrons without 
loosing gauge invariance. Unfortunately, there are problems and ambiguities in the the 
selection of form factors F 1>2 (see e.g., [53] for a detailed discussion). It has been shown 
already in Ref. [28] that NN bremsstrahlung contributions depend sensitively on the choice 
of these form factors and that cross sections calculated with no form factors are closest to 
the data. In our calculations, we have used the prescription [43,53] of using no form factors 
at the electromagnetic vertices of the nucleon term and replacing F 1 by the nucleon charge 
and F 2 by nucleon anomalous magnetic moment. 



D. Resonance- nucleon-photon vertex 



TABLE III. Coupling constants for the resonance-nucleon-photon vertices used in the calcula- 
tions. For iV*(1520) proton couplings are given in the first line and neutron couplings in the second 
line. 



Resonance 


9i 


92 


93 


A(1232) 


5.416 


6.612 


7.0 


iV*(1520) 


3.449 


5.074 


1.0 




-0.307 


1.862 


3.0 
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For spin- 1 resonance-nucleon-7 vertices, the form of the vertex functions is 



1 RN-y 



— 1- 



2mjv 



3i(£: 2 )7a + 



2m 



N 



-P\ + 



9s(k 2 ) 

2,171 jv 



{-k u g<* x + k^g vX )Z 



(19) 



where Z is 75 for the even parity resonance and unity for the odd parity one. Values of the 
coupling constants g±, #2, and #3 used in our calculations are shown in Table III. The first 
two are taken from [43], while the last one from [54]. The vertex function given by Eq. (19) 
fulfills gauge invariance [55] by the way of its construction. 

It should be noted that the vector spinor vertices [Eqs. (9)-(12) and (19)] should in 
addition be contracted by an off-shell projector Q av (z) = g av — |(1 + 2z) r y a 'y v where z 
is the off-shell parameter [43,56]. This operator describes the off-shell admixture of spin- 
\ fields [57]. The choice of the off-shell parameter z is arbitrary and it is treated as a 
free parameter to be determined by fitting the data. However, recently, the authors of 
Ref. [58] have proposed a different nNA interaction which leads to amplitudes where spin-| 
components of the Rarita-Schwinger propagator (see in section E) drop out, thus making 
the off-shell parameters redundant (see, e.g., [53,58] for further details). The full implication 
of this prescription on observables calculated within the effective Lagrangian model will be 
investigated in future. In the present study we work with the couplings given by Eqs. (9)-(12) 
and (19). 



E. Propagators 

In the calculation of the amplitudes, the propagators for various mesons and nucleon 
resonances are required. For pion, p meson and axial-vector mesons, they are given by 



G7(?) 



(q 2 



—1 



—1 



mi 



q li q v /q 2 



q 2 — m 2 



9' 



q 2 — m 2 A/ 



(20) 
(21) 
(22) 



In Eq. (22), the mass of the axial meson is taken to be very large, as the corresponding 
amplitude is that of the contact term. The propagators for uj and a mesons are similar to 
those given by Eqs. (21) and (20), respectively. 

The propagators for the nucleon and the spin-| resonance are 



G N (Pi) 



i{t>i+m N ) 
p 2 - m% 



(23) 



x 



i{t>i+ m R) 



Pt 



[m R -i{Y R /2)Y 

2 PM + 



3 ' 1 



3m 2 R 



3m 2 R 



(pfr-P-r) 



(24) 
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In Eqs. (24), Tr is the total width of the resonance which is introduced in the denominator 
term (p 2 — m 2 R ) to account for the fact that the resonances are not the stable particles; they 
have a finite life time for the decay into various channels. T# is a function of the center of 
mass momentum of the decay channel, and it is taken to be the sum of the corresponding 
widths for its decays to pion and rho channels (the other decay channels are considered only 
implicitly by adding their branching ratios to that of the pion channel): 

r_R = r^jv-Tr + r R _, Np (25) 

The partial decay width T R ^ Np is calculated in the same way as in Ref . [33] . This method is 
also identical to that used in Ref. [29] ; the only difference is that while Ref. [33] uses a fully 
relativistic expression, that employed in [29] is its non-relativistic counterpart. The partial 
width T R ^ N7T is taken as [29] 



"- nM - r ° [k^j ) { 0.25 + KW 1 ' (26) 



with 



kM = [(^ 2 - (m n + m N f)(p 2 - (m n - m^) 2 )] 1 / 2 

In Eq. (26), £ and T are 2 and 0.095 GeV for the A* (1520) resonance, and 1 and 0.120 GeV 
for A isobar, respectively. 

It should be noted that no width is included in the p propagator [Eq. (21)] as it corre- 
sponds to a exchange process between two nucleons where q 2 is in the space-like region and 
hence p meson has a zero width. However, this width is finite in the p propagator in, e.g., 
AA — > NNp — ► NNe + e~ process where the corresponding momentum is time-like. 



F. Amplitudes and cross sections 

After having established the effective Lagrangians, coupling constants and form of the 
propagators, we can now proceed to calculate the amplitudes for various diagrams associated 
with the AA — > AAe + e~ reaction. These amplitudes can be written by following the well 
known Feynman rules [50] and calculated numerically. It should be stressed here that the 
signs of the various amplitudes are fixed, by those of the effective Lagrangians, coupling 
constants and the propagators as described above. These signs are not allowed to change 
anywhere in the calculations. 

The general formula for the invariant cross section of the A + A — > A + A + e + e~ 
reaction is written as 

^ = , h, m K: ^ i - p ')'- 4/ '' 2 n ^. 

2 yj [(Pi • P2) 2 ~ m%\ l Z7r J o=i 

where represents the total amplitude, Pi and Pf the sum of all the momenta in the 
initial and final states, respectively, and p a the momenta of the particles in the final state. 
In Eq. (28) m e represents the mass of the electron. The term | \ already includes a sum 
over final spin and average over initial spin degrees of freedom of all particles. More details 
of the evaluation of Eq. (28) are given in appendix A. 
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FIG. 2. Calculated invariant mass distributions for dileptons produced in proton-proton (left 
panel) and proton-neutron collisions (right panel) at the beam energies of 1.04 GeV, 2.09 GeV 
and 4.88 GeV. The contributions of the NN bremsstrahlung (non-resonance), and the delta isobar 
and TV* (1520) resonance states are shown by dashed-dotted, dashed, and dotted lines, respectively. 
Their coherent sum is shown by solid lines. 



III. RESULTS AND DISCUSSIONS 

In Fig. 2, we show the invariant mass spectra for dilepton production in both pp and pn 
collisions at bombarding energies of 1.04 GeV, 2.09 GeV and 4.88 GeV. It can be noted that 
in all the cases the dominant contribution arrises from the intermediate states consisting 
of the A isobar resonance. In fact the total yields are almost equal to the contributions of 
the A amplitude alone. The pn cross sections are about a factor of 2-3 larger than those 
for the pp reactions even at the higher beam energy of 4.88 GeV. For beam energies of 1.04 
GeV and 2.09 GeV, the contributions of the A*(1520) terms are similar to those of the NN 
bremsstrahlung processes for the pp collisions, while for the pn case the latter are larger 
than the former. However, at 4.88 GeV, the N* channel prevails over the bremsstrahlung 
one in both cases. Any way, the A isobar production terms are still dominant even at this 
high energy. 
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FIG. 3. Calculated invariant mass distributions for dileptons produced in proton-proton (left 
panel) and proton-neutron collisions (right panel) at the beam energies of 1.04 GeV, 2.09 GeV 
and 4.88 GeV. Solid and dashed lines show the cross sections obtained by coherent and incoherent 
summations of the amplitudes of various processes shown in Fig. 2, respectively. 

The contributions of the nucleon bremsstrahlung process (graphs involving only intermediate 
nucleon lines) are considerably larger in case of pn reaction as compared to those for the 
pp case for all the beam energies. We do not observe the reversal of this trend at 4.88 
GeV as seen in the soft-photon model calculations of Ref [27] (where NN bremsstrahlung 
is larger for pp reaction as compared to that for pn one). On the other hand, if the off-shell 
behavior of the effective NN interaction is described by using a T matrix which includes 
the A degrees of freedom, instead of the meson exchange mechanism as employed in our 
model, the pp bremsstrahlung contributions turn out to be larger than those of ours at the 
bombarding energy of 1 GeV. Yet due to lack of calculations for the pn case, it may not 
be possible to make comments about the relative pp and pn bremsstrahlung contributions 
within this approach. In any case, the coherent sum of the A and bremsstrahlung terms as 
obtained within our model are similar to those calculated with the T matrix method. 

The role of the interference effects between various terms is investigated explicitly in 
Fig. 3 where we compare the cross sections obtained by coherent and incoherent summations 
of the amplitudes corresponding to various processes shown in Fig. 2. The former are the 
same as those shown by solid lines in Fig. 2. It is already noted in Fig. 2 that the delta 
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FIG. 4. Total dilepton production cross section in proton-proton (top) and proton-neutron (bot- 
tom) collisions as a function of beam energy. The nucleon, delta isobar and iV*(1520) contributions 
are shown by crosses, filled triangles and solid squares, respectively. The total cross sections ob- 
tained by coherent and incoherent summations of the corresponding amplitudes are shown by open 
circles and plus signs, respectively. 

contributions are slightly larger than the total ones at the larger mass end of the spectrum 
for all the beam energies in both the cases. In Fig. 3, one can further see that the interference 
effects are noticeable towards the larger mass ends of the spectra for beam energies below 
4.88 GeV in case of pp collisions. For the 4.88 GeV case these effects show up also at lower 
values of the invariant mass. On the other hand, they are relatively smaller for the pn case 
everywhere. 

Contributions of various terms to the excitation function (integrated dilepton cross sec- 
tion) are shown in Fig. 4 for the pp (top) and pn (bottom) collisions. The total cross section 
at 1.04 GeV is slightly larger than that at 1.27 GeV in both the cases. However, after this it 
rises monotonously with the beam energy. The delta isobar terms dominate the total pro- 
duction yields for all the beam energies. For collisions, NN bremsstrahlung and iV*(1520) 
terms are of similar magnitude except at the beam energy of 4.88 GeV where the latter term 
is larger than the former. However, in the pn case, the NN bremsstrahlung contributions 
are larger than those of the iV*(1520) at lower beam energies while the two terms contribute 
similarly at 4.88 GeV. Furthermore, the difference in the interference effects of 
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FIG. 5. Comparison of the effective Lagrangian model results for the invariant mass spectra of 
dileptons produced in proton-proton collisions, folded with the experimental filter, with the data 
of the DLS collaboration [37] at the beam energies of 1.04 GeV, 1.27 GeV, 1.61 GeV, 1.85 GeV, 
2.09 GeV and 4.88 GeV. Various curves have the same meaning as in Fig. 2. The solid circles show 
the experimental data. 

various terms in two case, should also be noted. Cross sections obtained by coherent (shown 
by open circles) and incoherent (shown by plus signs) summations of the amplitudes differ 
very slightly for all the beam energies in case of the pn collisions and at energies smaller than 
4.88 GeV for the pp case. This can be understood from the fact that the regions where the 
interference effects are visible in Fig. 3, contribute very little to the integrated cross sections 
for these cases. However, the interference effects do show up even in the total dilepton yields 
for pp collisions at 4.88 GeV beam energy. 

The cross sections shown in Figs. 2 and 3 can not be compared data of the DLS collab- 
oration. They have to be folded with appropriate experimental filter, which is a function 
of the invariant mass (M), transverse momentum (p?) and the rapidity in the laboratory 
frame (yiab) of the produced dileptons. We have used the DLS acceptance filter (version 
4.1) in the folding procedure. In Fig. 5, we present a comparison of the folded invariant 
mass spectra (which also include the final mass resolution [29]) and the data from the DLS 
collaboration. We note that at the beam energies of 1.04 GeV and 1.27 GeV the effective 
Lagrangian model calculations are able to describe the data reasonably well in the region of 
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FIG. 6. The calculated dilepton invariant mass spectra for proton-proton collisions in the beam 
energy range of 1.04-4.88 GeV as comparison to the DLS data [37]. The effective Lagrangian 
model results are shown by the dotted lines (the same as those shown by full lines in Fig. 4. The 
contributions of it , i] Dalitz decay and p° and oj direct decay processes are shown by dashed, 
dashed-dotted, open circles, and open squares, respectively. The incoherent sum of these cross 
sections with those of the effective Lagrangian model is shown by the solid line. Solid circles 
represent the experimental data. 

M > 0.1 GeV. However, at beam energies higher than 1.27 GeV, our calculations fail to 
reproduce the data, being lower than them by factors ranging from 2-20. Clearly, with in- 
creasing beam energy, other dilepton sources become important; these include Dalitz decays 
of hadrons (ir°, i] and a;) [23,29,30], direct decay of oj and p° meson to dileptons, processes 
leading to multi-hadronic final states [27] and two pion annihilation [60]. Nevertheless, it 
is encouraging that the present effective Lagrangian model is able to account for the dilep- 
ton yields in pp collisions at beam energies around 1.0 GeV. Therefore, this theory could 
be useful in making predictions for the dilepton spectra to be measured by the HADES 
spectrometer at GSI, for beam energies in this region. 

In order to study the role of other dilepton production processes not considered in the 
present effective Lagrangian model, we have included the contributions of the Dalitz decays: 
7T° — > 7e + e~ and r\ — > je + e~, and the direct decays p° — > e + e~ and oj — > e + e~. The details 
of these calculations are described in Ref. [29] and important relevant formulas are listed 
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in appendix B. The cross sections of these processes were incoherently summed to those of 
the total mass differential cross sections (solid lines in Fig. 5) of the ELM. The results are 
shown in Fig. 6 together with the DLS data. At the beam energies of 1.04 GeV and 1.27 
GeV, we have shown explicitly only the contributions of tc° and ir° and rj decay processes, 
respectively together with the ELM results as cross sections for the other decay channels 
are insignificant here. At incident energies of 1.61 GeV and 1.85 GeV we show in addition 
the contributions of p° direct decay also; the u direct decay is insignificant at these energies. 
However, at 2.09 GeV and 4.88 GeV, contributions of all the decay channels as mentioned 
above are shown explicitely. 

It is clear that now the data can be described reasonably well for all the beam energies. 
This comparison suggests that the lowest points in the mass distributions of the dilepton 
stem from the 7r° Dalitz decay. The 77 Dalitz decay is important at the intermediate masses 
for beam energies between 1.61 GeV to 4.88 GeV, while p° and u direct decay processes are 
important at higher mass ends of the spectra for these beam energies. In addition to these 
processes, the multi-hadron final state bremsstrahlung mechanism could also contribute 
significantly in the low mass region at the beam energy of 4.88 GeV [27]. We would like 
to mention that although several processes are included in our results shown in Fig. 6, the 
double counting problem is unlikely to be present there. The Dalitz decay and the direct 
decay processes of dilepton production are excluded from the effective Lagrangian treatment. 
Furthermore, we have not considered the processes like AA — > NR — > NNp° — > NNe + e~ 
(termed as the subthreshold resonance production). Therefore, problems of double counting 
in the vector meson direct decay contributions as suggested in Ref. [61], are not present in 
our calculations. 

IV. SUMMARY AND CONCLUSIONS 

We investigated the dilepton production in the nucleon-nucleon collisions at beam en- 
ergies in the range of 1-5 GeV within an effective Lagrangian model, which in proven to 
describe well the pion and kaon production in A A collisions. Most of the parameters of this 
model are fixed by fitting to the elastic NN T matrix; this restricts the freedom of varying 
the parameters of the model to provide a fit to the data. Along with the NN bremsstrahlung 
process, the model also includes the excitation, propagation and radiative decay of A(1232) 
and A* (1520) intermediate nucleon resonant states. The coupling constants at vertices in- 
volving resonances have been determined from the experimental branching ratios of their 
decay into various relevant channels. The interference terms among various amplitudes are 
included in the total T matrix. 

The reaction proceeds predominantly via excitation of the A intermediate state in the 
entire beam energy range considered in this work. The contributions of A* (1520) terms 
are relatively small. The AA bremsstrahlung contributions are also weak in comparison 
to those of the delta isobar except for the higher mass region of the invariant mass spectra 
at the beam energy of 1.04 GeV. In this region of these spectra considerable interference 
between various terms is also visible. For the case of pp collisions, the AA bremsstrahlung 
and A* (1520) contributions are similar in magnitude at lower beam energies. However, for 
the pn case, the former is larger than the latter at these energies. On the other hand, for 
the beam energies of 4.88 GeV, A* (1520) terms are larger than those of the bremsstrahlung 
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in both the cases. A key result of our study is that the pn bremsstrahlung is stronger than 
the pp one even at the beam energy of 4.88 GeV. 

At the lower beam energies (1.04 GeV and 1.27 GeV), the effective Lagrangian model 
calculations are able to explain the data of the DLS collaboration for the invariant mass 
distribution of dilepton emitted in pp collisions except for the lowest mass point. This is 
encouraging in the context of the analysis of the new experimental data on the dilepton 
production expected shortly from the HADES collaboration at GSI, Darmstadt. 

However, at higher beam energies (1.61 GeV, 1.85 GeV, 2.09 GeV and 4.88 GeV), the 
effective Lagrangian model underpredicts the DLS data. This indicates that with increasing 
energy, other dilepton production sources become important. To stress this point further, 
we incoherently summed the 7r° and rj Dalitz decay and p° and uo direct decay contributions 
to the cross sections of the ELM. With this procedure, it is possible to provide a good 
description of the DLS data on the invariant mass distribution at all the beam energies. This 
also removes the discrepancy observed between the effective Lagrangian model calculations 
and the DLS invariant mass distribution data for the lowest mass point at the beam energies 
of 1.04 GeV and 1.27 GeV. 

Our study has clarified the role of A and iV*(1520) intermediate baryonic resonance states 
in the dilepton production in the NN collisions. The interference effects of various terms 
are visible towards the large invariant mass ends of the spectra for proton-proton collisions. 
However, the DLS experimental filter effectively cuts off this mass region. Therefore, these 
effects do not show up when calculations are compared with the data. This may imply 
that an incoherent combination of resonance and decay contributions is sufficient to explain 
the DLS data. However, with improved quality of the data this impression could change. 
Therefore, precise data where also the hadrons in the final channel are detected would be 
welcome. This would put constraint on the models calculating the role of various mechanisms 
in the dilepton production in nucleon-nucleon collisions. The present effective Lagrangian 
approach should be improved by including the amplitudes for the decay processes also so 
that they can be summed coherently with the amplitudes of the other channels considered 
already. 
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APPENDIX A: SOME DETAILS OF THE PHASE SPACE AND THE 

AMPLITUDE 

In this appendix we give additional details of the calculation of the four-body phase 
space factor and the amplitude A* % . 

The kinematical situation for calculating the differential cross section [Eq. (28)] is illus- 
trated in Fig. 7 where pi and p 2 depict the initial momenta of two nucleons while p 3 and 
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P4 their final momenta. The momenta of e + and e~ are shown by p + and p_, respectively. 
P represents half of the total center of mass (CM) dilepton momentum. The differential 
cross section can be written as 



da 



4m%m 2 e P 2 Q 2 L Q 2 N 



dMdtt P dPdtt QL d{l QN [E l+ Q L - E^Pcosdr] [{E 3 + E A )Q N + (E 3 - E 4 )Pcos6 2 ] 

X E l+ E p \P p \(2n)* lAff ' (A1) 

where Ql = p_ - P and Qtv = P3 + P- &\ and 62 are angles between P and Q^, and P 
and Qat, respectively. In Eq. (Al), we have defined Ei + = E + + E_ and E\_ = E + — E_, 
where E + and E^ are energies of e + and e~, respectively. E p and \P P \ are the energy and 
momentum of the projectile nucleon in the CM system. The invariant mass of the dilepton 
is given by, M 2 = (p + + p ) 2 . Integrations over the residual degrees of freedom in the 
differential cross section [Eq. (Al] are carried out by using the Monte Carlo techniques. 




FIG. 7. Illustration of the kinematical situation for e + e~ pair production in NN collisions, pi 
and P2 are the initial nucleon momenta in the CM system while P3 and P4 are the same for the 
final nucleons. P is half of the total dilepton momentum. 

A total number of 90 diagrams [42 for NN bremsstrahlung, 16 involving A intermediate 
states and 32 involving iV*(1520)] contribute to the amplitude A- fl . Due to their large 
number, it is hard to calculate them by usual trace techniques. Instead we follow the 
procedure, given below. Similar method has also been used in our previous calculations 
within this model. 
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We first note that amplitude corresponding to each individual diagram can be split into 
a hadronic part and a leptonic part as 



Aft ~ tf^l* (A2) 

where k is the photon momentum. For example, for the diagram of type (a) of Fig. 1, 
involving pion exchange and the nucleon intermediate state (NN bremsstrahlung) , we have 

h f t = _ { (9nnt£\ 1 ^ s 4 )75^(P2, s 2 )u(p 3 , s 3 )rJ r7V7 
^ \2m N J q l -m 2 K * 

{i)i+m N ) 

x — 2"75gw(pi, Si)/ iso , (A3) 

Pi ~ m N 

lf % = eu(p + , s + )e-f u u(p_, s_), (A4) 

where Sj) are the free Dirac spinors with four momentum pi and spin Sj, and is the 
same as defined in Eq. (16). fi SO is the isospin factor which for each graph, is obtained from 
the separate treatment of the isospin parts. In Eq. (A4) e is |(1 + T3). Similar expressions 
for the amplitudes corresponding to all the graphs can be written in a straight forward way 
by using the Lagrangians and the propagators given in the main text. 

The summation over spins of the modulus square of the leptonic part gives 

T _ V /^t//» 

— £^ <- M <v 

spins 

e 4 

= ^2 [- 2 9^M 2 + A(p + ^ u + p +I/ p_ M )J (A5) 
In the evaluation of the hadronic part of the amplitudes we encounter terms of the form 

T? 1 = u( P3 , g3 ) 7 / ( ;^ + ^ rf / >(p 1 , Sl ), (A6) 

where is a linear combination of 16 4 x 4 Dirac matrices (1, 7^, cr^, 75, 7s7 At )and rnj 
is mass of either nucleon or the resonance. In the latter case it also includes the resonance 
width as shown in Eq. (24). Further calculations of such terms are carried out by using the 
manipulator package REDUCE. The matrix elements so calculated are summed together as 
complex numbers. The result still depends on the initial and final spins of the nucleons. The 
modulus square of this sum is then averaged over the initial spins and summed over the final 
spin states. This leads to |A^| 2 . By this procedure a large number of matrix manipulations 
are avoided. 



APPENDIX B: CALCULATIONS OF DALITZ DECAYS AND DIRECT DECAY 

OF MESONS 

In this appendix we describe briefly the calculations of Dalitz decays r\ — > 7e + e~ and 
7T° — > 7e + e~, and direct decays V — > e + e~ (V = p and u). More details can be found in 
Ref. [29]. 
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The calculation of the process pp — > MX — > / ye + e~ proceeds in two steps. First the 
meson M (rj and ir) is produced in the pp collisions and then the Dalitz decay M — > / -/e + e~ 
is considered. In the procedure of Ref. [29], the prametrizations of Refs. [3,14,62] are used 
to calculate the M production cross sections. The rj Dalitz decay to / ye + e~ is given by (see 
also [63]) 



j2\ 1/2 

x(l-^) |iW e -(M)| 2 , (Bl) 



dM 3n M \ M 2 J \ M 2 ) 

where the form factor is parameterized as 

F v _ ne+e -(M)=\l-—j . (B2) 

The value of parameter is 0.72 GeV. In Eq. (Bl), a represents the fine structure constant, 
and r^_ +2 7 is the partial width for the rj — > 27 decay whose value is adopted from the particle 
data group compilation [64]. 

The Dalitz decay 7r° — > r )e + e" is calculated in the same way using the following form 
factor 

F^ ie+e - (M) = (l + B n oM 2 ) , (B3) 

with B n o = 5.5 GeV^ 2 . 

Similarly, calculations of the direct decay of p and uj mesons proceed in two steps; the 
production of these meson in pp collisions is calculated in the first step and their decay to 
dilepton in the second 

da(M) NN ~> VNN ^ e+e ~ NN = da(M) NN - VNN r v ^ e+e -(M) 

dM = dM rY ot (M) ( } 

The calculation of the first term on the right hand side is described in details in Ref. [29] 
and will not be repeated here. The decay width T v ^ e+e - (M) is given by 

777 ^ 

rw +e -(M)=cv^, (B5) 

where Cy = 8.8xl0~ 6 and 0.767xl0~ 6 for V = p and uj, respectively. The total vector 
meson decay width Y^( ot is defined in the same way as in Ref. [29]. We refer to [61], for a 
more recent description of these calculations. 
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